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Abstract 

Aluminum-oxygen reaction is important in many highly energetic, high pres- 
sure generating systems. Recent experiments with nanostructured thermites 
suggest that oxidation of aluminum nanoparticles occurs in a few microsec- 
onds. Such rapid reaction cannot be explained by a conventional diffusion- 
based mechanism. We present a rapid oxidation model of a spherical alu- 
minum nanoparticle, using Cabrera- Mott moving boundary mechanism, and 
taking self-heating into account. In our model, electric potential solves 
the nonlinear Poisson equation. In contrast with the Coulomb potential, 
a " double-layer" type solution for the potential and self-heating leads to en- 
hanced oxidation rates. At maximal reaction temperature of 2000 °C, our 
model predicts overall oxidation time scale in microseconds range, in agree- 
ment with experimental evidence. 
* karen . mart irosyan@ut b . edu 



Many previous studies have been directed to understand the mechanism 
and kinetics of aluminum particle oxidation. [1] -[9]. Aluminum oxidation 
exhibits high enthalpy and has been extensively used for propulsion, py- 
rotechnics and explosion reactions [10] [H]. Nanoenergetic materials (NM) 
based on aluminum thermites may store two times more energy per vol- 
ume than conventional monomolecular energetic materials [12]. The size 
reduction of reactant powders such as aluminum from micro- to nano-size 
increases the reaction front propagation velocity in some systems by two to 
three orders of magnitude [12] [13]. The development of novel NM, their 
design, synthesis and fabrication procedures are critical for national secu- 
rity and it was recognized to be significant to support advanced weapons 
platforms. Among numerous thermodynamically feasible Metastable Inter- 
molecular Composites mixtures the most widely investigated are Al/Fe203, 
AI/M0O3, AI/WO3, Al/CuO, Al/BisOs and AI/I2O5 nano systems M'M- 
The main distinguishing features of these reactive systems are their signifi- 
cant enthalpy release and tunable rate of energy discharge, which gives rise 
to a wide range of combustion rates, energy release and ignition sensitivity. 

There are several advantages of using Al/Bi203 and AI/I2O5 nanocom- 
posites: (i) reduced ignition and reaction times; (ii) superior heat transfer 
rates; (iii) tunability of novel energetic fuel/propellants with desirable phys- 
ical properties; (iv) enhanced density impulse; (v) incorporating nanoener- 
getic materials into the MEMs and NEMs systems [12]. Our recent experi- 
ments suggest that oxidation of nanoparticles of aluminum with Bi203 and 
I2O5 occurs in a few microseconds [I7]-[19]. Rapid reaction in the nanostruc- 
tured thermites cannot be explained by a conventional mechanism based on 
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the diffusion of Al and O atoms in oxides. 

In tliis report, we present a rapid oxidation model of splierical aluminum 
nanoparticles surrounded by oxygen, using Cabrera-Mott oxidation model 
[7]-[9] witli a self-consistent potential, and taking self-heating into account. 
Using nonlinear self-consistent potential in Cabrera-Mott model, as opposed 
to the Coulomb potential, is more accurate, and yields higher oxidation rates. 
For nanosized particles, nonlinear model gives "double-layer" type solution 
for the potential, since potential changes rapidly near the metal-oxide inter- 
face. As a result, the electric field at the metal- oxide interface, determining 
the oxidation rate, is less sensitive to oxide thickness, and stays nearly con- 
stant at intermediate stages of oxidation (while for the Coulomb potential, 
electric field decreases, and the oxidation rate quickly drops). This nonlinear 
effect for the potential, combined with self-heating, leads to rapid tempera- 
ture increase, resulting in a dramatic (orders of magnitude) increase of oxi- 
dation rates for nano-sized particles. Our model can predict oxidation times 
in microseconds range, in a good agreement with experimental data [12]. 

To estimate the reaction times, we assume that aluminum sphere (radius 
25 nm, with thin oxide layer of 3 nm) is surrounded by oxygen. The sphere 
is rapidly heated to ignition temperature Tq, sufficient to initiate oxidation 
reaction, further boosted by self-heating as a result of oxidation. We assume 
the spherical symmetry of the problem. 

In the Cabrera-Mott model of metal oxidation [I]- [3], [7] aluminum ions 
are helped to escape aluminum boundary (overcoming ionization potential 
W) by a self-consistent electric potential V. The electric field is induced in 
the oxide layer by imbalance between excess positive aluminum ions and 
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electrons in the oxide layer, due to difference in chemical potentials of the 
metal and oxidizer. Excess concentration of electrons and ions in the oxide 
layer depends on the electric field potential V via appropriate Gibbs factors, 
leading to a self-consistent version of Poisson equation for V, which in the 
spherically-symmetric case is given by 

= ('''^^) " ^""^"^ ^^'""^ (^) ' ^1-^-^2, 

(1) 

Here ri is the metal particle radius, (r2 — ri) is oxide layer thickness, ko — 
8.99 lO'^ e = 1.60 IQ-^^C, N = {n^ni)^''^ , and and are concen- 

trations of respectively aluminum ions and of excess electrons in the oxide 
layer. 



rie = 2 (27rmeA;6r//i2)'/'exp 
( eW^A 



Ui — Ni exp 



where vrif. is the mass of electron, fcf, Boltzman constant, h Plank constant, 
Ni is the concentration of sites available for hopping metal ions. Thus 

7V = Kn.r = iVo(j)%xp(-^), 

A^o = (2Ar.)5 {2T^mekbTa/h^f'^ - 1.5 10^^ m-^ at Tq ~ 7bQK. 

The physical meaning of Wi is the difference of chemical potentials for metal 
ions in the metal and the oxide; is the potential difference for electrons in 
the conduction bands of aluminum metal and the oxide (a semi-conductor); 



and the value of Vq is determined from the condition that metal ion con- 
centration at the interface with the oxide equals rij. In Cabrera-Mott model, 
those ionization potentials may be considered as model parameters [Z]-[9]. 

Oxidation of metal occurs via tunneling of aluminum ions into the ox- 
ide layer, overcoming ionization potential of maximum height W > Wi, and 
assisted by the self-consistent electric potential. Electric field provides poten- 
tial energy decrease for ion hopping from the bottom to the top of ionization 
potential, at a distance of ao ~ 0.4 nm. This leads to an equation for the 
metal-oxide interface radius ri : the normal velocity of the metal boundary 
Un is determined by the electric field {—W) at the boundary, appearing in 
a Gibbs factor. 



Here rio is the initial metal sphere radius, which we assume to be 22 nm, 
Qi ~ 0.0166 nm'^ is the volume of oxide per aluminum ion, n2 ~ 10 nm~^ is 
the number of metal ions per unit surface area, u ~ lO^^s"^ is frequency of 
tunneling attempts, g = 3 is aluminum valency. 

In the Cabrera-Mott model, it is assumed that escaped metal ions migrate 
to the outer boundary of the oxide where they react with the oxygen, while 
local Gibbs distribution of excess densities of electrons and metal ions inside 
the oxide is essentially unaffected. Due to the spherical symmetry, the radius 
of the oxide-oxidizer interface r2 changes uniformly, and can be found from 
conservation of the number of metal ions, taking into account difference in 





ri(0) = rio. 



5 



volumes per metal ion in the metal and the oxide, (r| — t^q) = —n {r-\ — rfg) , 
K ^ 0.386. Thus 

1 

rs = rs (ri) = [rl^ + k (r?o - r?)) 3 , ^ ^ 0.386. (4) 

Here r2o = '"lo + ^i where 5 ~ 3 nm is the initial oxide layer thickness. 

For small metal particles, it is important to take self-heating into ac- 
count, due to heat released by exothermic aluminum oxidation, resulting in 
temperature increase of the remaining metal and oxide layer. For nano-sized 
particles, temperature can be assumed to be uniform. Thus temperature can 
be computed based on reaction heat release and specific heats of reagents. 
Assuming constant specific heats, 

T^T(n)=To + P^.{rl,-rJ)^^ 

CaiPai^I + Cai^03Pa1203 V 2 '^IJ 

Here p^;, PAI2O3 are densities of the aluminum and oxide, c^,, c^i^oa 
specifics heats per unit mass, H^i is the oxidation reaction enthalpy per unit 
mass of aluminum, and a is the proportion of released heat which is used up 
for self-heating. We take reaction initiation temperature Tq = Ib^K. Since 
r2 can be found from ri using Q, the temperature T in ([s]) is determined by 
ri. Due to high enthalpy release in aluminum oxidation (~ 24 kJ per gram 
of aluminum ), adiabatic assumption a = 1 yields unrealistically high maxi- 
mum temperature, even when melting and vaporization heats are taken into 
account. We assume that a = 0.11 of the heat released contributes to self 
heating, while the rest is lost due to radiation, heat conduction, and convec- 
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tion. For such value of a, the maximum reaction temperature, corresponding 
to ri = 0,r2 = r2(0) in Q, Q, will be Tm ~ 2000"C, which agrees with 
experiments [T2] . 

We use experimentally determined maximum temperature to control the 
effect of the heat transfer process on oxidation rates. Detailed modeling of 
the heat transfer is outside the scope of this Letter. However, we note that 
an assumption of fairly large heat loss is sensible. For a nanoscale particle, 
the convection loss mechanism is important. Indeed, in the steady state 
limit, heat loss of a small sphere due to convection can be modeled by the 
Newton's law of cooling, with the heat transfer coefficient r/, [Wm^'^K^^], 
given hj 1] = ^ [21j, where fc, [Wm~^K~^], is air's thermal conductivity, 
and r is Al sphere radius (this limit is applicable for small Grashof numbers, 
which is the case for nano-spheres) . Thus, the convective loss rate ri{ATTr'^)6T 
is proportional to r, while the total heat released due to chemical reaction is 
proportional to r^. As a result, convective heat loss becomes significant for 
small spheres. An estimate shows that for an aluminum particle of radius 
25 nm, all the heat released by oxidation may be lost to steady convection 
in a time scale of 10~* s, which is comparable to the time scale of the final 
rapid phase of the oxidation process in our model. Radiation loss, which 
is proportional to r^T'^, contributes less than convection for the values of 
parameters in our model. 

We note that Eq. (|3| for the metal boundary velocity is similar to the 
model of metal sphere oxidation considered in [IJ , however we use a different 
potential (a solution of the nonlinear Poisson equation ([T]) rather than the 
Coulomb potential), and we also take the self-heating effect into account. As 
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a result, oxidation rates are dramatically increased. 

We have found the self-consistent potential V, solving numerically the 
boundary value problem for the Poisson equation ([T| for various radii ri (we 
have used the Newton's method to solve a discrete version of the problem). 
Our results for the potential are illustrated in Figure 1. Our values of ion- 
ization parameters Vq = 0.65 V, (p + Wi = 1.5 V, W = 1.7 V, a = 0.4 nm are 
consistent with the data in the Cabrera-Mott paper [7j. We used rio = 22 
nm for the initial metal radius, and = 25 nm for the initial metal+oxide 
radius. When the remaining metal radius is ri < rio, the radius of the metal 
-|-oxide r2 can be found from Eq. (|4]), and the temperature, with the self 
heating taken into account, is determined by (|5]). We use Tq = 750K as the 
initial temperature. 

We note that for the case of cylindrical symmetry, the potential equation 
is equivalent to Painleve 3 (with particular values of parameters). 

Solution of the nonlinear Poisson equation ([l| enables us to compute the 
gradient of the potential on the metal surface, and corresponding potential 
decrease due to electric field near metal-oxide interface, for various radii ri. 
In Figure 2, we show electric potential decrease at a distance a = OAnm 
from the metal-oxide interface for the solution of nonlinear Poisson equa- 
tion (solid line), and compare it with a potential decrease computed using 
Coulomb potential (broken line). Potential decrease determines oxidation 
speed, via an exponential Gibbs factor in Eq. ([3]). Our results demonstrate 
that the nonlinear model for the potential, with the self-heating effect, yields 
considerably higher oxidation rates. 

Using results for the gradient of the potential on metal-oxide boundary, 
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we can find the radius of tlie metal spliere as a function of time, by solving 
the moving boundary Eq. ([s]). This equation is separable, since it follows 
from Eqs. (|4]),([5]) that the right hand side is a function of ri only. Thus, a 
solution can be found by a numerical integration. The solution is shown in 
Figure 3. The overall oxidation time scale is in microseconds, much faster 
than for macroscopic particles. We further observe that in our model most 
of oxidation occurs very quickly towards the end of oxidation process, since 
reaction rate dramatically increases with the temperature rise due to the 
self-heating. 

Our model yields oxidation time of 32 fis for a spherical aluminum nanopar- 
ticle of radius 25 nm at initial temperature of 750 K. By using exactly the 
same ionization parameters and the self-heating mechanism, but taking the 
Coulomb potential rather than a solution of the self-consistent equation ([T]), 
oxidation time is calculated to be 2.5 ms, that is almost 100 times slower. 
As for the diffusion limit, it is known to correspond to the linear in the 
Coulomb potential V approximation of the Eq. ([3]). It is clear from Eq. ^ 
that the oxidation rate in the linearized model will much slower than for the 
Mott model. Using the same model parameters and the self heating as in 
our model, oxidation time in the diffusion limit would be about 40 ms; and 
without the self-heating it is estimated to be 14 seconds at T = 750i^. Exper- 
imentally, an explosive oxidation reaction of nanosized particles is initiated 
at a temperature of about 750i^ (i.e below aluminum melting temperature), 
with the reaction time scale estimated to be in microseconds [12]. This 
comparison strongly suggests the Cabrera-Mott model with a self-consistent 
potential as the most likely reaction mechanism. 
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Figure 4 illustrates the oxidation rate of change of mass of aluminum in 
a metal-oxide sphere with initial radius 25 nm , as a function of aluminum 
mass and temperature T. Such a rate of change of reacting mass appears e.g. 
in Semenov explosion model, with usually given by the Arrhenius factor 
exp(— ^) with a constant E. However in Cabrera-Mott oxidation model, E 
depends on the remaining mass and temperature (since electric potential 
drop in the self-consistent Cabrera-Mott model depends on the metal radius 
and temperature). 

Nonlinearity in equation for the self-consistent electric potential, together 
with the self- heating effect, lead to significant increase of oxidation rate for 
nano-sized aluminum particles, compared to the model with Coulomb po- 
tential and without self-heating. Results of our modeling suggests oxida- 
tion time scales in microseconds range for nano-sized aluminum particles, in 
agreement with the experimental evidence [T2]. Results of computation of 
oxidation time scale are very sensitive to values of ionization potentials, since 
they appear as Gibbs factors in the equation for boundary velocity. There- 
fore we expect that oxidation rate will be very sensitive to defects, and to 
deviation from spherical symmetry. 

After this work was completed, we have learned that in a recent paper 
[1] a different mechanism of oxidation initiation, based on high temperature 
run-off and the Coulomb potential, was proposed for a smaller, 10 nm radius, 
particle with very thin initial oxide layer. In our model with the nonlinear 
equation for the potential, rapid full oxidation of nanosized particle can occur 
at lower, experiment-matched, maximal temperature, due to nonlinear effects 
for the potential near metal-oxidizer interface. 
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Figure 1: Self-consistent potential V for ri = 2, 6, 10, 14, 18, 22 nm. Here rio = 22 
nm, rao = 25 nm, Vq = 0.65 Y, (p + Wi = 1.5 \ , W = 1.7 V, Tq = 750 K. 
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Figure 2: Potential drop a distance a = 0.4 nm on the oxide side from aluminum 
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metal-oxide interface, shown as a proportion of the total Cabrera-Mott potential 
Vq = 0.65 V, as a function of the remaining metal radius, and taking self-heating 
into account. Solid red curve: nonlinear self-consistent model for the potential; 
black dashed curve: Coulomb potential. Here rio = 22 nm, r2o = 25 nm, 0-1- Wj = 
1.5V ,W = 1.7 V, To = 750 K , Tm = 2273 K. 
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Figure 3: Oxidation time scales: (a) ratio of remaining aluminum metal mass to 
initial aluminum mass, as a function of time. Here rio = 22 nm, r2o = 25 nm, 
(f) + Wi = 1.5 Y , W = 1.7 V, To = 750 K , Tm = 2273 K, a = 0.4 nm. (b) final 
stages of oxidation shown; solid curve: ratio of the remaining aluminum metal 
mass to the initial aluminum mass, dashed curve: ratio of temperature increase to 
the maximal temperature increase 
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Figure 4: Oxidation rate of change of mass of metal aluminum in a metal-oxide 
sphere of initial radius 25 nm, as a function of metal mass m and temperature 
T (K), in logj^o scale. Metal mass is shown as a ration to the initial mass of the 
particle mo. Here initial radius is 25 nm, (j) + Wi = 1.5 V , W = 1.7 V, To = 750 
K, Tm = 2273 K, a = 0.4 nm. 
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